Superfluid ^He in globally isotropic random media 
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Recent theoretical and experimental studies of superfluid ^He in aerogels with a global anisotropy 
created, e.g., by an external stress, have definitely shown that the A- like phase with an equal spin 
pairing (ESP) in such aerogel samples is in the ABM (or, axial) pairing state. In this paper, 
the A-like phase of superfluid ^He in globally isotropic aerogel is studied in details by assuming a 
weakly disordered system in which singular topological defects are absent. Through calculation of 
the free energy, a disordered ABM state is found to be the best candidate of the pairing state of 
the globally isotropic A-like phase. Further, it is found through a one-loop renormalization group 
calculation that the coreless continuous vortices (or, vortex-skyrmions) are irrelevant to the long- 
distance behavior of disorder-induced textures, and that the superfluidity is maintained in spite of 
lack of the conventional superfluid long range order. Therefore, the globally isotropic A-like phase 
at weak disorder is, like in the case with a globally stretched anisotropy, a glass phase with the 
ABM pairing and showing superfluidity. 
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I. INTRODUCTION 

Superfluid ^He in aerogel, which is a random medium 
with a high porosity, has been studied as a prototype 
of impurity scattering effects on an anisotropic Cooper- 
paired system [ij. The aerogel has a structural correla- 
tion length ^a, corresponding to a typical distance be- 
tween neighboring strands, of the order of 30 - 50 (nm) 
which is comparable with the pairing coherence length 
^0 — fivp I (2'KkBTc) in the pressure range relevant to the 
superfluid '^He. For this reason, the scattering events of 
the quasiparticles due to the aerogel structure are differ- 
ent from those of electrons in amorphous-like dirty metals 
corresponding to the situation with ^ Co and seem to 
be characterized by a local anisotropy of the scattering 
amplitude [l|,[3,[i. 

In recent years, the presence of a high pressure phase 
with an equal spin pairing (ESP), called the A-like phase, 
has become an active research subject, because it had 
been believed previously that the B-like phase with the 
Balian-Wherthamer (BW) pairing is the only stable su- 
perfluid phase in aerogel However, several NMR 
experiments have suggested the presence of the A-like 
phase near Tc{P) and a strange lowering of the polycriti- 
cal pressure (PCP) accompanying this phase 5, 6]. Even 
theoretically, there have been some suggestions favoring 
the presence of the Anderson-Brinkman-Morel (ABM) 
pairing state at least at short scales |7| by assuming a 
local anisotropy of aerogel structures, although the low- 
ering of PCP has not been discussed there. On the other 
hand, it has been argued that the ABM pairing state in 
aerogel has no superfluid response at long distances, as 
a consequence of random orientations of 1- vector pinned 
by the local anisotropy of the aerogel and the resulting 
proliferation of nonsingular coreless vortices, or vortex- 
skyrmions [1, where the unit vector 1 is the orbital 
anisotropy axis and expresses the direction along which 
the energy gap vanishes. However, experiments seem to 
show nonvanishing and anisotropic superfluid responses, 
like in the bulk liquid Another pairing state, called 



the robust pairing state, was proposed as a candidate 
of the A-like phase showing superfluidity. However, 
it was difflcult to identify this state, which is not ther- 
modynamically stable in the bulk liquid, as the A-like 
phase. 

This controversy on the pairing state of the A-like 
phase has been resolved in more recent studies 0, 
14] for the cases where the aerogel has a global anisotropy 
brought by an external stress. NMR measurements in 
both of uniaxially compressed [l^l and stretched 
aerogels have been nicely explained by assuming that 
the A-like phase in these aerogels is in the ABM pair- 
ing state with a proper alignment of the orientation of 
1-vector. On the other hand, it seems that there is no 
consensus at present on the pairing state of the A-like 
phase in the globally isotropic case [1, H, [lH . Even if the 
A-like phase in this case is also in the ABM pairing state, 
the fundamental question on the presence or absence of 
superfluidity in this case 0] needs to be resolved. 

In the present work, thermodynamic stability of the 
ABM pairing state in globally isotropic random media is 
examined in details through comparison on the free en- 
ergy between different pairing states, including the pla- 
nar and robust states, by assuming some of real aerogels 
to be globally isotropic. Further, the presence of the 
quasi long range superfluid order in such a disordered 
ABM state is established at weak disorder where singu- 
lar topological defects are absent. Here, the superfluid 
order will be reasonably defined through the correlation 
function [8] 



G(R) = ReTr(Ap(r + R) A;(r)) 



(1) 



between the spin-triplet gap parameters ■^p(r) and 
Ap(r + R). Here, the zero temperature limit is assumed 
so that the thermal fluctuation of the gap parameter may 
be neglected. Further, the overbar denotes the random 
average, and Tr expresses both of the trace in spin space 
and the average over the relative momentum p on the 
Fermi surface. This gap parameter, which is a tensor 
in spin space, depends not only on the ordinary ampli- 
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tude and phase but also on the orientations of spin and 
orbital degrees of freedom of Cooper- pair s. At larger 
scales than the dipole coherence length [l^ , the spin ori- 
entation is locked in the orbital one corresponding to the 
I's orientation, and a short range correlation of the I's 
orientation corresponds to a short range superfluid order 
measured by the correlation function ([T]). The presence 
of a quasi long range superfluid order suggests that the 
corresponding superfluid correlation length is inflnite. At 
a glance, one might wonder that such a long range corre- 
lation be destroyed by vortex-skyrmions which are gener- 
ated by continuous textures of the 1-vcctor. However, we 
find based on a renormalization group (RG) analysis that 
the vortex-skyrmions appearing at short scales in glob- 
ally isotropic systems may be irrelevant perturbations at 
long distances, implying that a nonvanishing superfluid 
response is well-defined. Therefore, the A-like phase at 
weak disorder is expected to show superfluidity, just as 
seen experimentally [l^l- A brief sketch on the free en- 
ergy calculation in the present work has been reported 
elsewhere [1] previously. 

In sec. II, the Ginzburg-Landau model including effects 
of randomness is derived in a form useful for a free energy 
calculation, and the free energy is evaluated in details in 
sec. Ill based on the Gaussian variational method often 
used in random systems. In sec. IV, the presence of a 
quasi long range superfluid order is explained by per- 
forming one-loop diagram calculations accompanying a 
functional RG method, and results are summarized and 
discussed in sec.V. Some of technical or numerical details 
will be explained in Appendices. 



II. DERIVATION OF GINZBURG-LANDAU 
ACTION IN DISORDERED CASE 

As a starting microscopic Hamiltonian for deriving a 
Ginzburg-Landau action or functional, we choose the 
BCS Hamiltonian with an attractive interaction in the 
purely p-wave channel, which is written in the familiar 
notation as 



h = kn = I unit, where 



Hp - fiN = ^ £,p aj,^„ flp^Q + Hint, 



(2) 



where 



^i„t = -3\g\Y,Ol,iq)d^,,i<i), 



OMj (q) = X] T^^-P+q/2,a(i '^2)0/3 ap+q/2,/3-(3) 
p '^PP 



Performing the standard decoupling [la] in Hi^t by in- 
troducing the pair-fleld A^j^ j, where n (j) denotes the 
3-components of the spin (orbital) degree of freedom, 
the superfluid part of the partition function is given by 

(T,exp[-/„^/^dsi7int(s)]) = / VAVA* exp{-S) in the 



n 



2 ^ 

q 

h.c, 



(At(q))^,„ 



ds a. 



'p+q/2 



a{s)a. 



■p+q/2,l3{s) 

(5) 



/p denotes the momentum integral J(Pp/{27T)^, 

(Ap(q))a_^ = Af_,_i{q)pi{iaf_,a2)a,(3 is the pair-fleld, 
i/'cr(r) = X]p ap,o-e'P '' is the quasiparticle field, and ( ) ex- 
presses the ensemble average over the quasiparticle dis- 
tribution. The GL action is obtained by, in S, keep- 
ing just the quadratic and quartic terms in A^^j. The 
pair-field is assumed to be independent of the imagi- 
nary time s because quantum fluctuations of A^.; do 
not have to be included in considering superfluid phases 
of ^He in equilibrium, in which fluctuation effects are 
safely negligible. For '^He in aerogel, the total quasiparti- 
cle Hamiltonian needs to include a term associated with 
an impurity scattering. As usual, it will be expressed 
hereafter as a nonmagnetic random potential term 
Wimp = J2a I d^i' u{r)ipl{r)'ipcr{r) . The scattering poten- 
tial u(r) has zero mean, and the quasiparticle life time r is 
deflned by the relation t^^ = 27rA^(0)(|up_p'p)p, where 
A^(0) is the density of states per spin in the normal state, 
and ( )p denotes the angle-average over the orientation 
of the relative momentum p on the Fermi surface. If the 
aerogel we assume has no global anisotropy, t deflned 
above is independent of the quasiparticle momentum p'. 
Using a quasiparticle Green's function [3] G£(p, p') de- 
fined prior to the impurity average, the quadratic part 
S2 of S is expressed as 



S2 = T 



-^E 

q,q' 



^E 



pjp' 



PiPj 



X cap + q/2, p' + q'/2) G-e(-p + q/2, -p' + q'/2) 



X A*,(q)A^,,(q'). 



(6) 



In the present situation where the critical fiuctuation 
is negligible, a q-dependence of A^^- follows from the 
quenched disorder. In the GL regime where the am- 
plitude of A^,j is small, it is sufficient to keep, in 52, 
disorder-induced terms related to the q-dependences of 
Afj^ j, and the corresponding contributions from the GL 
quartic term may be neglected. Then, the quartic term 
in our GL action takes the same form as the familiar one 
for the disorder- free bulk liquid '^He (see, e.g., Ref . (T5| ) 

54 = (.Pl\A,,^A,^^\^ + P2{A;^.A,,^y 

+ P^Ap^ ^A^ ^Afj^jA^j + /34A^ jAjy^iA^j-Apj- 
+ /35A;,,A,,,A;^.A,,,). (7) 
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(b) 



FIG. 1: Diagrams expressing Gor'kov boxes leading to eqs.(9) 
and (10). The solid line with arrow and the dashed line de- 
note the quasiparticle Green's function and the impurity line 
carrying r~^, respectively. 



In the weak coupling limit without any vertex correction, 
54 is obtained from the expression 



Si 



^ / (Gap)G_,(-p))2Tr(AtA^A^Ap)(8) 



expressing the Gor'kov box of Fig. 1(a), where G^{p) = 
(is — ^p)~^ with e = e + sgn(e)/(2r) is the impurity- 
averaged quasiparticle propagator. The resulting 54. wc is 
given by replacing (3j in eq.(I7]) by (3^^'^\ where (^^2"'^^ = 



-PI 



(wc) 



2/3p^ = 2/3(*'=)(T), 



N{0) 



E 



2407r2T2 ^ (2n + 1 + (27rrT)-i)3 

n— \ / / 



E 



7C(3)„^„(2n 



l + (27rrr)-i)3' 



(9) 



and C(3) ~ 1.2. 

In performing the impurity average in eq.(l8]), the lead- 
ing term in (Epr)^^ was kept by neglecting diagrams 
with crossed impurity lines |18| . We need to comment on 
our neglect in eq.([H]) of two types of vertex corrections 
induced by the impurities. First, the impurity-ladders 
dressing the four corners of Fig. 1(a) were neglected. 
These vertex corrections of the order of 1/ {2ttTt) are 
present even in the p-wave pairing case because we take 
account of wave vector-dependences of the squared im- 
purity potentials litkp. However, they merely renor- 
malize the magnitude of (3o and never affects a rela- 
tive difference between free energies of two different pair- 
ing states. On the other hand, irrespective of the pair- 
ing interaction, there are also additional diagrams, de- 
scribed in Fig. 1(b), accompanied by a single impurity 
line. These diagrams do not contribute to /jj™'^'' and 



/3 



while they change other coefficients as follows 



0,[T1 : p: 



(wc) 



-/3o^imp, P. 



>c) 



and Pc^ 
{2ttTt) 



(wc 



p, 

1«1. 



(wc) 



Pt^ 



-PoA 



imp J 



/3oAiinp, where Aimp is of the order 



As is well known, the coefficients Pj (j = 1 to 5) ap- 
pear in a manner dependent on the pairing state in the 
condensation energy of bulk ^He. Hereafter, such Pj^s 
combinations in the ABM, BW, planar, and the robust 
states are denoted by Pabm, PbWj Pp, and /3r, respec- 
tively 20] which will be indicated as /3n in the lump 
hereafter. In the disordered case, they are expressed as 

Pabm — P245 — /3o^imp, 
Pbw — Pi+ P2 + — Po\mp, 

Pp = Pi + P2 + ^ - PoAi^p, 



pK = /32+ g(/3i3 + 5/345) -/3oAi,„p, 



(10) 



where Pi- 



P, + Pj, and P^jk = A + Pj + Pk- Note 



that Aimp appears in the same form in all /3ns. Further, 
a pairing state with a small /3n has a lower free energy 
at the mean field level. Therefore, relative stability be- 
tween the different pairing states cannot be reversed by 
including the contribution of Fig. 1 (b) . 

Besides, Pj must include the so-called strong couphng 
(SC) corrections [21,22,19] which, in clean limit, sta- 
blize the ABM state as the bulk A-phase of superfluid 
■^He. In the disordered case, the SC corrections to Pj 
consist of two contributions. One is the expression in 
clean limit with the Matsubara frequency e replaced by e. 
The other consists of terms including impurity-induced 
vertex corrections. Hereafter, they will be denoted by 



Sp^-"^'^ and 5Pf^' , respectively [19|. That is, we have 



□ (sc) 



P, = P 



(wc) 



5P 



(sc) 



5~pf^\ 



Details of Spf""^ and (5/3]'"^ 



were examined in Ref. 19] thoroughly. Their pressure de- 
pendences in each pairing state are needed in obtaining 
a theoretical phase diagram, and their numerical values 
will be illustrated for reference in Appendix A. 

Now, let us turn to detailing the sec- 
ond term of ^2 by expanding fij{£) = 
S^Sp'P^P'pe{v+.v'+)G-e{-V-,-p'_) in powers of 
the impurity potential u, where p± = p ± q/2. Up to 
0(u2), we express as flf + /f i^-' + flf , where 



G,(p+)G_,(-p_ 



(11) 



and 



- / p,p,G,(p+)(G_,(-p_)G_e(-p+ + q') 

G-,(-p-)Ge(P-+q'))^^q-q', (12) 
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p Jk 



G,(p + k/2)Ge(p - k/2)G_,(-p - k/2)G_e(-p + k/2)uk+qU,k-q. 



p + - j 1 Ge(p + k/2)(G_,(-p - k/2))^G_,(-p + k/2) 



+ (G^P + k/2))^G,(p - k/2)G_e(-p - k/2) ]uu+qU^u-q' 

r 



(13) 



The contributions in ^2 corresponding to eqs.® and 
(jlip give the ordinary quadratic term in the so-caUed 
Abrikosov-Gor'kov approximation and in the weak- 
coupHng hmit 0, . Its expression is well known and 
given by 



(0) 



where 



-'1,3 



- ( 2Kiq^qj + (5jj 



(14) 



N{0) 
3 

K2 = 



ln( ) + V^(l/2 + l/[4^rT]) - ^(1/2) 



(15) 



Tco is the superfluid transition temperature of the bulk 
liquid, ')jj{z) is the digamma function, and 



^0 



2t:T\ 12 



^E(n + l/2 + l/(4^rr))- 



(16) 



is the coherence length. 

In fij, the first order term J^e^iP^^) easily found 
to vanish after performing the p-integral. Thus, we have 



only to focus on /, 



(2) 



(2') 

After substituting /.^' into eq.®. 



a larger |k| is found to become dominant in the result- 
ing replicated action Sdis (see below), while, for |q|, 
|q'| < 2tt$,q^, any q and q' dependences included in the 



p-integral are negligible compared to the large |k|. Then, 

(2) 

the p-integral in f^j is derived in the conventional man- 
ner [Tsj used for obtaining the static supefluid response, 
and we obtain 



k-q' 



(17) 



Although, by substituting this into 52, a disorder- 
induced term appears in the GL action, it is more conve- 
nient to directly work in the corresponding quartic term 
iSdis arising after the impurity-average of the free energy, 
where the index "dis" imphes "disorder". To do this, let 
us introduce the replica-description of the averaged free 
energy Q 



F = -T hm 



1 





(18) 







FIG. 2; Diagrams giving Sdis 



where = Trexp(— 5). The quartic term 5dis, de- 
scribed in Fig. 2, appears in the replicated GL action S 
and is given by [25 1 



a,fc=l \ i- I I \ // a,b quq2,q3 

^ ^ ^ ^ ^ a,fc 91,92,93 

X A(,'^].(g2)AW(9l+93"g2), (19) 



r 



where '0'^^^(z) — dip{z)/dz, and It is easy to verify that 5dis can also be obtained by 



dk 



l^irN{0)y I ^^l"k|4. (20) 
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assuming the following quadratic action to be present in 
the original action S : 

52,dis = y"d'V[C/(r)5,j +y(r)(5,,, -3a,(r)a,(r))] 

X A;jr)A^,,{r) (21) 

Here, yields a Gaussian ensemble satisfying = and 
3ai{r)aj{r) — 5i j, while the potentials U and V have zero 
mean and satisfy U{r)V{r') — 0, and 

U{r)U{r') = V{rjV{?)^T-'^5l3dS^^'>{r-r') (22) 

with 

**-*Wj|i5(«''''(i + 3;;i^))^ (23) 

In this way, one can regard the original GL action S 
below Tr as 



S — ^o"^ + ^4 + 52,dii 



(24) 



Hereafter, the pair-field ^ will be expressed by sep- 
arating the amplitude |A| from the symmetry variables 
[m composed of the spin and orbital degrees of freedom 
together with the overall phase 4> in the manner 



(25) 



Following the standard notation, a^^i in the ABM state 
is given by 



(m + in), 
V2 



(26) 



with the triad (m, n, 1) of unit vectors. On the other 
hand, it takes the form, e'* Rf_i.i/^/3 and e'* Rf_i,k ^fk/^^ 
for the BW and planar states, respecti vely , where Rf^^i is 
a rotation matrix, and 6jj = 5i,j — |15j . According to 
the definition of the 1-vector mentioned in sec. I, the same 
notation on the anisotropy axis will be used for both of 
the ABM and planar states. Then, 5dis is rewritten for 
the ABM pairing state in the form 5dis(i) +'5dis(2) i where 

5dis(i) = -^l^Y. j '^'r-|A('^)(r)A('')(r)p, (27) 



a.b 



and 



dis(2) 



1 3 

2rio 



(5/3d^ / dV|A('^)(r)A(^)(r)|^ 



X [(l('^).l(''))2_l]. 



(28) 



The corresponding action for the planar state takes the 
same form as above, while S^is is given, in the BW and 
robust states, simply by 5iSdis(i)/6. Since eq. (|^5|) is non- 
vanishing only if the 1-vector is spatially varying so that 
(1(a) . \{b)f _ 1 is nonzero, |A(")(r)A(^) (r)|2 in eq.i 



may be replaced by its mean field value |Amf|^ below 
the critical region in the close vicinity of Tc, or as far 
as a slowly varying a^_j is assumed. Of course, |AmfP 
needs to be determined by examining F . 

Before ending this section, it will be appropriate to 
discuss about the treatment on the impurity scattering 
used in this paper. Our procedure on the impurity scat- 
terings used in and ^4 is more or less an extension of 
the Abrikosov-Gor'kov approach [2^ based on the Born 
approximation. Except in the situation with extremely 
weak disorder, the so-called unitary limit including mul- 
tiple scattering processes is often used by assuming the 
isotropic s-wave scattering event to be dominant (see, 
e.g., Ref.iT,])- On the other hand, an inhomogenuity of 
the order par ameter amplitude |A| to be created spon- 
taneously [26[ by impurity scatterings was not incorpo- 
rated in the present analysis. In an isotropic approxima- 
tion, this effect becoming more important at higher pres- 
sures was studied in an elaborate numerical work [27| . 
Throughout the present paper, however, we argue that 
the local or global anisotropy of scattering events in aero- 
gel, which has not been incorporated in calculations in 
previous microscopic works [7|, |27|, is indispensable for 
describing the features in the phase diagram associated 
with the A-like phase of liquid '^He in aerogel. Further, 
to examine effects of the local anisotropy, one needs to 
derive an expression of a disorder-induced term, corre- 
sponding to 52,disj in the GL action. In order to achieve 
these purposes consistently, we have chosen to work in 
the simplest Born approximation and its extention. To 
perform a more quantitative comparison between exper- 
imental and theoretical phase diagrams and obtain re- 
sults on physical quantities such as the temperature de- 
pendence of |A| comparable with experimental data, the 
multiple scattering events and spatial variations of |A| 
need to be incorporated within a model of anisotropic 
and random scattering. 



III. FREE ENERGY AND GRADIENT TERMS 

To evaluate free energy for various pairing states in the 
present disordered case, we will use the Gaussian varia- 
tional method (GVM) . In this method, a trial Gaussian 
ansatz Sg for the replicated action S is first invoked. 
Then, the total free energy F is evaluated as 



F = F„ 



-iS-S 
n 



(29) 



where Fg is the free energy for Sg divided by the number 
of replicas n, ( )g is the ensemble average on Sg, and 
the n ^ limit is taken at the end. The GVM has 
been satisfactorily applied in evaluating free energy of 
the random Ising-spin 28] and elastic systems [2^ . 

To apply GVM to the present problem, we will first 
examine how to determine an appropriate trial action in 
our case with a couple of fields, the amplitude fluctuation 
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6\A\ = |A| — |Amf| and a^j- consisting of the symmetry 
variables. Since we are not interested in a negligibly nar- 
row critical region in the close vicinity of Tc, we will as- 
sume, as usual, the two variables S\A\ and a^j- to be sep- 
arable in the trial Gaussian action. This assumption on 
the trial action greatly simplifies our analysis for eg. ([29]) . 
In fact, the Gaussian approximation does not have to be 
assumed in the original action S which appears only as 
its average in eg. ((29)) . To clarify this point, let us rewrite 
the original gradient term 



2,grad 



included in S2^\ It is not difficult to see that the Ki 



term in eq. (|30p is rewritten as 



2T 



^1 {dial^idja^,^ + dja*^^idiaf,^j) 



Re{a;^,a,,,){2d.\A\d,\A\-d,d,\A\') 



(31) 



except surface terms. In the ABM or planar state, the 
presence of a^^i in the 2nd line of ea. (PT|) makes this term 
a nonGaussian form, because the factor Re(a* ja^j) be- 
comes {dij — lilj)/2 there, although it is merely 5ij/3 
in the BW or robust state from the outset. In the dis- 
ordered ABM or planar state, however, the 1-vector has 
no orientational long range order Q, and hence, the ran- 
dom average of the factor Re(a* jO^j) is merely Si_j/3 
irrespective of the correlation range of the 1-orientation. 
In this way, the original gradient term of S, if applied to 
ea. ((29)) . can be replaced by 



5. 



2,grad 



K{V\A\f 



\A 



(32) 



for all pairing states considered in this paper, where 
K = {3K2 + 2Ki)/6. Here, according to the assump- 
tion of a slowly varying j mentioned below ea. (j28p . 
the factor |Ap was replaced by its uniform value |AmfP 
to be determined later. 

Then, in the total and averaged free energy F — 
-Famp + -f^sym, the (5|A|-part Famp and the a^^ - part 
Fsym can be treated independently below: 



yn 1 

Ti T amp 
i'amp = -J hm . 

n— +0 n 

2 

^'sym = lim "^'^ 



Tts\a\ exp(-5amp), 

Trj;,, , exp(-5 



(33) 



(34) 



where A^^^i = |AMF|a/i,i- Since variations of af_i,i are al- 
ways accompanied by Amf in '52,grad, the free energy 
correction due to the purely thermal fluctuation of sym- 
metry variables is independent [30] of |Amf| and thus, of 
the details of pairing states. Since such a free energy cor- 
rection insensitive to |Amf| should take a common value 
to all p-wave pairing states , this purely thermal correc- 
tion will not be considered in {S)g hereafter in examining 
a relative stability between different pairing states. 

According to treatments performed so far, the repli- 
cated action iSamp for the ABM state is given by 



n ^ 

^amp = T ^ y ^ / d' 
a=l •' 



lA^-^^P + i^(V|A(" 



)n2 



+ /?ABM|A(«)|4-^5/3d^|A(-)|2|AW| 



6=1 



(35) 



The corresponding expression for the BW (robust) state 
is given by replacing /3abm and the factor 3/5 in the sec- 
ond line by /3bw {Pk) and 1/2, respectively, while the 
corresponding one in the planar state follows from re- 
placing /^ABM by ^p. _ 

On the other hand, the replicated action, 5sym = 
^giad + 5^113(2): for the ABM and planar states is 



S 



\A 



grad 



MF| 



2T 



E 



2KidiauSja 



dis(2) 



-I- K2dia^jdia*^^ 



^<5/3d|AMF|^E((l^''^-l^''^)'-lX36) 



fc=i 



if the field a^j in the planar state is represented by 
eq.(Cl) in Appendix C. 

Here, for later convenience, the gradient energy in the 
purely ABM pairing state will be expressed in the hydro- 
dynamic representation [Tsl . [sij 



■ad 



_ If" 



26Cv(") • L^'^) + 2Cv('') • curll^'^) 



(37) 



where 



Vi = mjV.rij, L = l(l-curll). My = Si^j-Ahlj, (38) 

with positive constants A and 5, and Spr is the replicated 
Frank energy term 



1 r " 



(39) 
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in the terminology of the nematic hquid crystal, if the 
1-vector is identified with the nematic director. 

On the other hand, for the the BW and robust states, 
the (5 /3d term of eg . ([5^ is absent. Since, as mentioned 
earlier, the thermal fluctuation term of the symmetry 
variables is unnecessary for the present purpose, even 
{Sgiad)g in eg. ipg)) may be neglected. Therefore, for the 
BW and robust states, we have no contribution of -Fsym, 
and the total free energy F can be identified with -Famp- 
Now, let us turn to evaluating free energy of the dis- 
ordered ABM state. The corresponding results for other 
pairing states will be commented on at the end of this sec- 
tion. First, to examine i^amp, it is convenient to rewrite 



eq. ()35|l in the form expressed in terms of a scalar order 
parameter 0(r) 



i(a)^(fc))2 



(40) 



which was studied within GVM in Ref.f28| as a contin- 
uum model of a random Ising spin system. Here, the 
scale transformation, | A|2(A")3/2/(r(A^(0))i/2) ^ ^2/2 
and [iV(0)/i<r]^/^r r, was performed. Details of deriva- 
tion of the free energy for the model ((40|) is explained in 
Appendix B. By rewriting eq. (|B12p . the resulting Famp 
is found to take the form 

Fa„.p (^(0))\2 T(7V(0))3/2p„ ^ , 



V 



4/3abm " 
— ■ — \iq — lu)\ — 



27ri^3/2 27r' 

Vc x2 



where V is the volume. 



a 



N{0) 



(3g 



(41) 



(42) 



and 



TP. 



ABM 



(Af(0))l/2K3/2' 

3<5/3d 



5/3, 



ABM 



-9- 



(43) 



The dimensionless momentum cutoff Pc/(2tt) will be as- 
sumed below to be unity. We note that the third term 
of eq.(|41|) merely gives a negligibly small correction to 
the flrst and second ones in the relative difference be- 
tween free energies of two different pairing states, since 
g and u are at most 0{T^/Ep). Depending on the disor- 
der strength, this correction may be negligible compared 
with the contribution from Fgym which will be examined 
below (Note that Fgym is absent in the BW and robust 
states) . 

In contrast to Fampj it is not tractable to directly eval- 
uate Fsym in the ABM state. To evaluate Fgym in a 



different manner, let us first start from examining free 
energy of the simpler model [2^ 

+ r-2^^(l-cos[2(e'(") -6l(''))])]. (44) 

b 

Assuming a Gaussian trial action 

= ^EE^aV(q)f^''^Hq)f?''H-q), (45) 

q a,b 

the corresponding averaged free energy Fxy is given by 



FxY 

TV 



1 

2n 
W 



c I q'TTg{q)+ / Trln[g-i(q)] 

q 



VV X — ^ 

— (^exp(-2i?,,(0)) + n) 



(46) 



except a constant term, where Bab{0) — /q(5aa(q) + 

Gbb{<l) — 25ah(q)), and the n ^ hmit is taken at the 
end. By following the procedures in Ref.^2^, the disorder 
dependent term of the first term of ea.(P5|) is given by 



^ du [(7]u 



(47) 



while the integrand in its second term is expressed by 



Trln [^-^(q)] = n 



;;„2n 



In (eg 



du ( \a 
— m' 



(48) 



Details of the function \a\u can be seen in Ref.[29'|. Using 
the properties of [(t]„ carrying the disorder strength W ^ 
the q-integral of the second term of eq. can be shown 
to be convergent. Then, it is found that, up to the lowest 
order in the disorder strength, the sum of the first and 
second terms in ea. (H51) is disorder- independent. There- 
fore, the change of free energy density induced by the 
quenched disorder is, up to 0(V7), given by the second 
line of eq. ((46|) . i.e.. 



FxY(VK)-fxY(0) 

V 



w 

2T 



exp 



1 



(49) 



which is independent of T because c is an elastic con- 
stant divided by T. Note that eq. (j49p is determined by 
the behavior at short scales of 0{p~^), implying that the 
free energy is unaffected by the details of long distance 
behaviors [32|, i.e., the presence or absence of quasi LRO. 
In fact, refiecting the fact [1^ that the elastic behavior at 
short scales is determined within the replica-symmetric 
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approximation, the result p9|) coincides with the corre- 
sponding one of the random-force model [H, [s^l 



Srf = / d\ 



-(Vef + f{r)9{r) 



where f — 0, and 



fir)f{r')=AT-^Wd''^\i 



(50) 



(51) 



This action is equivalent to the Gaussian replicated ac- 
tion obtained from eq. (|44p with the replacement 1 — 
cos(2(6i('^) - 6*^)) ^ 2(6'('^) - 6i(''))2. 

Based on this fact for the model (|44|) . we have evalu- 
ated -Fgym by, in the last term of Sgymi keeping only the 
lowest (harmonic) order terms in Euler angles represent- 
ing the 1-vector. Then, if using the representation 



I = z cos9i + (x COS0/ 4- y sin(/)i) sin6'/, 



(52) 



one finds that the resulting last term of iSgym takes the 
form 



>dis(2) 



2or 



,5/3d|AMF|'^ 



g(f')^2 



r (53) 



which depends only on 9^°''^ and flp"*. Further, the gradi- 
ent energy, ea. ([57)) . in the ABM state will be replaced, 
for simplicity, by its isotropized version 



-.(iso) 



erad 



1 

2r 



Hp 

a=l 

_l_ 2C(''^°) v(°) • curll(°) 



(iso)j^^(a)^2 



s 



Fr 



(54) 



corresponding to the limit of a Bose gas of molecules 
with the ABM pairing symmetry [sHi, where and 

^(iso) 

are averaged coefficients which follow by replacing, 
e.g., 1 • vl • (curll) in the original action by (ljj)ui(curll)j 
and applying the absence of the 1-orientational LRO. The 
second term, proportional to v • curll, will be neglected 
hereafter because it simply becomes a sum of purely sur- 
face terms after expressing it via the Euler angles. This 
easily follows from the fact that, in the representation 
([5^ . V • curll is proportional to (Vcos(26';) x W4>i)z — 
2(Vcos^0; X Vsin0;)a; -|- 2(Vcos^0; X Vcos0;)y. Further, 
the remaining terms expressed in terms of the Euler an- 
gles, 9i and 0/, will also be linearlized by using the ab- 
sence of LRO. For instance, sin(26';)(V0i)^ will be re- 
placed by (sin(26';))(V6'/)^ which vanishes due to the ab- 
sence of LRO. The resulting expression is Gaussian in 
V0; and V6'/, and there are no cross terms like V6'/V(/); 
there. In fact, the Gaussian term in V0;, which is the 
relevant one for the present purpose (see eq. fSS]) ). results 
only from Sfi-. In this way, the relevant gradient energy 



term in Ssym becomes 



a—1 



(55) 



in the weak coupling approximation where Kh = SKg = 
3Kt = 3|AMFp(i^i+i^2)/4. The coefficients ^C,,, K^, and 
Kf including the strong coupling corrections are given, 
up to the lowest order in (Tc — T)/Tc, by their weak 
coupling expressions divided by the mass enhancement 



factor [36], if Sf3j^'^' and Sp'^^'^'' are properly incorporated 
in l3j appearing in |Amf|- Thus, eg. ([55]) is expected to 
be applicable even at higher pressures as far as pressure 
dependences of and A''(0) are incorporated through 
their experimental data. The remaining (/);-dependent 
terms are purely thermal fluctuation contributions unre- 
lated to the quenched disorder and hence, may be ne- 
glected hereafter to derive the (5/3d-dependent correction 
to the free energy. Then, eq. ([55)1 accompanied by eq. ([55]) 
is of the same form as the random- force model, eg. ([50)1 . 



if 3T^/3d|AMF|^/20 is identified with W, and hence, the 
resulting disorder contribution to -Fsym is given by 



sym 



('5/3d)-i^sym(0) 



V 



~9TN{0)\Xp\Sl3i 
257r/3ABM(A'i + if2)^o 2^' 



(56) 

We are now at the stage of discussing stability of the 
pairing states and the resulting phase diagram of super- 
fluid ^He in globally isotropic aerogels. To perform the 
remaining task, we need just the free energy expressions, 
egs. ipTj) and ([SS]) . and information on the SC effects in 
each state (see Appendix A and Ref.p^l. First, judging 
from the data of SC parameters, there is no possibil- 
ity that the ABM state is replaced by the robust phase 
|ll| . The contributions from the SPd term to the free en- 
ergy definitely show that this term favors the anisotropic 
ABM and planar states. Although the disorder effect on 
the SC parameters may suggest a small gain of the con- 
densation energy in the robust state, it is quite difficult 
to, in the weak disorder regime, find such a situation that 
the robust state is realized due to an enhanced disorder. 
Rather, it is more reasonable to examine the planar state 
as a candidate, other than the ABM one, of the A-like 
phase. However, since inevitably f3p > Pabm, ^amp in 
the planar state is higher than that of the ABM state. 
In addition to this, the planar state is not favored even 
through -Fsym: As shown in Appendix C, the gradient en- 
ergy in iSsym of the planar state is 2.4 times bigger than 
that of the ABM case. Since the expression for the planar 
state corresponding to eg. ([5S|) is also inversely propor- 
tional to the magnitude of the gradient energy, the free 
energy gain in the planar state due to the random sym- 
metry variables is much smaller than that of the ABM 
state. By taking account of these results on F altogether, 
we conclude that even the planar pairing state cannot be- 
come stable as the A-like phase in the GL region in ■^He 
in aerogels. 

In Fig. 3, a typical example of pressure v.s. tempera- 
ture phase diagrams we obtain is shown. There, egs. (|4ip 
and ((56|) were used based on data of pressure dependences 
of and Tcq. The disorder-induced reduction of PCP 
indicated as a solid circle is a consequence of the large 7- 
value used here. The parameter 7 was defined in eg. (|20p 
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Temperature (mK) 



FIG. 3: Example of calculated P-T phase diagrams obtained 
based on the present free energy analysis. The solid curves 
denote the transition curves in the disordered case specified 
by the parameter values {2nTTco)~^ = 0.058 for P = 20(bar) 
and "/ — 21, while the dashed ones are the corresponding bulk 
transition curves. The hatched region indicates the A-like 
phase region. The Tco{P) and Ef{P) data are taken from 
Table 4.1 of Ref.fH] and Table VI of Ref.[33|, respectively, 
while the analysis on the SC correction entirely follows the 
phenomenological method in Ref . jT9| . 

as a measure of the strength of the local anisotropy in the 
scattering events. The fact that the A-like phase is lim- 
ited to such a narrow temperature range is a combined 
effect of Tc-reduction and the SC effect ^] shrinking 
with increasing disorder. 

IV. QUASI LONG RANGE ORIENTATIONAL 
ORDER IN A-LIKE PHASE 

In the preceding section, a typical phase diagram fol- 
lowing from evaluation of free energy was shown in Fig. 3. 
However, it is important to note that, at the present 
stage, the transition curves in the figure merely imply 
changes of the pairing states. As noted in sec. I, if the A- 
like phase of ^He in aerogel is in a disordered ABM pair- 
ing state, the genuine long range order of 1-orientation 
is absent in the A-like phase Q. It is crucial to clar- 
ify whether such a 3D phase with no genuine long range 



superfluid order may show superfluidity (see sec. I). In 
this section, we address this possibility at weak disor- 
der where the singular topological defects are not excited 
via the disorder. This issue is highly nontrivial because, 
in the globally isotropic case, the nonsingular vortices 
[38| or vortex-skyrmions [1| may appear as a result of 
a disorder-induced 1-texture and thus, may destroy su- 
perfluidity. It is shown below that a one-loop renormal- 
ization of the gradient e nergy terms accompanying the 
functional RG treatment [39ll40l |4l| yields only a stable 
fixed point at which the vortex-skyrmions are irrelevant. 
This implies that the A-like phase at weak disorder may 
have quasi long range superfluid order and superfluidity. 

To examine long distance behaviors of the symmetry 
variables of the disordered ABM pairing state, we exam- 
ine i^sym again. For the sake of the ensuing analysis, 
however, the kinetic part of 5giad expressed in the form 
of eq. (|37p will be rewritten in the form 

Po{l~A) J 

where denotes / (Pr, and v = v — 6CL/[po(l — 
A)]. Next, eliminating the longitudinal component of 
V, V in eq. (|57p is replaced by its transverse component 
/ d^r'{V X (V X v'))/(47r|r - r'|), and eq.^ can be 
replaced by 

5g..ad - 5f,. ^^£^JJ^ J^ X (V X ^^^Hr))h 

X Mlf{ri)w{r - ri)w(r' - ri)(V' x (V x v(")(r')))j 

where w{r) — |r|^^, and v' denotes v(r'). Further, using 
V^w(r) — ~4:7rS^^\r) and rewriting the terms in eq. ([58|) 
consisting only of V X L^"^ we obtain the following ac- 
tion of the nonlocal gradient energy appropriate for the 
ensuing RG analysis 



pn^^^r) ■ 0(")(r') - y [fl^^^r) • (V x L('^)(r')) + ft'^'^^r') ■ (V x L^'^^r))] 



^[n'^'^^r) ■ (V X (V X + f2(")(r') • (V x (V x l^"''))] - ^divL(")(r) div'L('^)(r') 



E 



•"1 a=l 



2f 



[(7(f2('^)(r) X V),(f2(")(r') x V)^ - ai [((V x L(")(r)) x V),(f2(") (r') x V'), 



+ ((V' X L('^)(r')) X V)^{n^''\r) x V)^] -I- era divL^-^) (r)V, div'L^'^) (r')V;.] u;(r - ri)w(r' - n). 



(59) 
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where 

n,(r) = (V X v(r)), = Eyfel • (V,l x Vfcl) (60) 

is the Mermin-Ho relation [38] in the ABM state, Sijk is 
the antisymmetric tensor, and the following redefinition 
of the Frank energy term 

1 r " 

= ^ /E(^Ml^''^^A.l^"'+^2(divlW)2 
^ a— 1 

+ A3((l('^) • V)l('^))2) (61) 

has been done. Further, the relations cJ^l^S^l^ = 
(divl)2 + + ((1 • V)l)^ and P = 1 were used. Note 
that the coefficient Kt of the twist deformation term was 
absorbed into T to define T. The term, which is absent 
in the bare action, has been included because it is gen- 
erated via renormalization. The bare values of the seven 
coefficients except Jii in eg. ([59)1 are positive, although 
their detailed expressions are not necessary in our analy- 
sis given below. In fact, it will be assumed that, through 
the dipole energy term, the d- vector c?^ (see ea. ipS)) ) is 
locked into 1 at large scales of interest so that the gradi- 
ent term on may be absorbed into the Frank energy. 
Nevertheless, we have the stability conditions 

As + 1 > 0, A3 + 1 > 0. (62) 

The goal in this section is to find an action at a stable 
disorder-induced fixed point by examining the scaling of 
the coefficients. 

Following Ref.js^ in which a functional RG analysis 
was performed for ^pr, let us focus on T limit, in 
which thermal fiuctuation effects arising from higher or- 
der terms in T are neglected, and determine the form 
of 5giad at the stable fixed point. To perform this, the 



disorder energy term will be generalized to 




In the functional RG analysis based on e — 4 — 
d-expansion, a stable disorder-induced fixed point is 
determined by R{z) of 0(e) in magnitude and the 
fixed point values of the coefficients in 5grad- To 
perform the one- loop renormalization of 5giad, the 
1- vector with the momentum q of 0(1) in mag- 
nitude will be expressed in terms of the "trans- 
verse" variables (pj in the momentum-shell (e~' < 

|q| < 1) as l(r) ^ T(r)y^l-E,=i,2('^'^'Hr))2 + 

Ej=i.2 0^^'nr)e(^'Hr), where e(j') • 1 = e^^) • e^^) = 0, and 
the disorder function R{z) will be expanded in powers 
of 0*^^^. Further, when examining a renormalized iSgrad 
up to one loop order, the replica-index dependences of 1 
and e'-') may be neglected |39| . Then, we only have to 
examine the correction — {dSgrad SSdis) to 5giad in T 
limit, where (^iSgiad is the second order correction in c/)^^^ 
to 5grad, and 

1 r " 
<55dis = -i^ /^^'Hl)E E((^^^'^)^'^nf)r, (64) 

where R^^\l) — dR{z) / dz\z=i, and a term which van- 
ishes in n ^ limit was neglected [s^ . 

To illustrate the one-loop renormalization procedure, 
let us first focus on the isotropic limit with A = h = Q 
in which the original gradient energy is given by ea. (j54p . 
Alternatively, one may start from S-pr + {p/T) J^, w{r — 

r')Ea^^"Hr)-^^^"''(r') in place ofeq.(P|. For simplicity, 
the replica index a and its summation will be omitted 
hereafter. By noting that the v • cur 11 term has no bulk 
contribution and thus, is negligible, we find 



SS. 



■(iso) 
grad 



w{r - r') --^((0O))2 + ((0(^))')2)O . n' - 4^^(^-)</.('=)e«a,Jpef (T x da)xs,.M'^ 

+ 4^(9^0(^)1 • (e(^) X aj) - d,<p^^n ■ (e(j") x 8^1)) a;(0('=))'l' • ((ef'^))' x dj') 

+ ^1 Et^'p '^'^Ip^'^SplA - 94.94.5,- + A2[(T • V)Ip(I • V)lAe(^)ef - (divT)^^,,,] 

3.k 

+ A3[(e(-'") • V)I,(e('=) • V)I, + e«(T • V)Tpef (T • V)Ta - 2((T • V)!)^^,- fe] ]<p^i^^^\ 



(65) 



r 



where the remaining terms harmonic in 
1 



A5] 



Fr 



2T 



^[<5,.fc[(V0W)2-|-A3((T-V)0 



can be identified with the " noninteracting" action for the 
(/.(^■)-fields. _ 

I evaluating — {5S gra-d^S dm) ^ we encounter the foUow- 



j.k 



+ A2((e(^) • V)</)(-'')(e('=' • V)0(*^)] 



(66) 
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ing expressions in the momentum-shell 



/(r) 



E 



((0«(f))20O)(r),/)('=)(r))A,fe(r), 



dence of the "noninteracting" action, eq. ((66)) . on these 
coefficients. In all terms in the one-loop renormalization, 
however, the result of integration in the momentum-shell 
is expressed by the quantity J(A2,A3), and its depen- 
„ „ dence on A2 and A3 is found not to affect the result- 

/ / ^(r — r')c)p9^((0'*^(f))^(/)'-'^(r)0'^'^^(r')) ing fixed points and the hnear stabihty around them. 

JrJr' Thus, tn simnlifv the ensin'np' exnressinns_ the denen- 

dence 



to simplify the ensuing expressions, the depen- 
B (r- r' — r) (67) and A3 will be omitted hereafter. 



After the trivial integration in the momentum-shell, we 



T^E^A^^{r)J{X2,X3)il-e-'), where 
1 — e~' is the thickness of the momentum- shell. Here, 
the A2 and A3 dependences of J arise from the depen- 



easily obtain / 



Then, using V w(r) = - 

4^rVo<5^,.E.5-(r;0)(i 



47r^'^'(r), we find /p,,y(r) = 
e-')/3, where Jq = J(0,0). 



Therefore, using the relations J^j ^, 
U) ~ ^\e^fdJx]l M, we have 



and dfj,e 



\er'd„lx]l 



{SSgradSSd 



-(l-e~') 



,,i?W(l)Jo 



T 



2p 



w{r~r')n{r) ■ n{r') - IGtt 1 



[(1 - A3)a4.a^L + 2A2(divT)2 + (4A3 - A2)((T • v)T)2] 



(68) 



Taking account of the rescaling factor e'*^'^ of T p3l |. 
we obtain the following recursion equations 

^f-i = f-i(2-e- Joi?«(l)(l-A3-/5)), 

^ = -Joi?^i)(l)A2(l + A3 + /5), 
dl 

^ = -Joi?(i)(l)(3A3-A2 + A3(A3 +/})), 
dl 

^ = -Joi?(i)(l)p(l + A3+p), (69) 

where p = 32ttp(1 — 1/d). The first equation simply 
ensures that, within the present analysis, we stay in the 
parameter space at zero temperature with no thermal 
fluctuation. Under the stability condition l-f A3 > 0, the 
following two fixed points 

(i) a; = \l = p* = 0, 

(ii) a; = Xy2, /}* = -1 - A* < (70) 



are found. The case (i) expresses the nematic fixed 
point [3^ with no vortex-skyrmions which is easily shown 
through a linear stability analysis to be a stable fixed 
point. On the other hand, the case (ii) expressing a 
fixed line has a negative value of p. However, this neg- 
ative value does not imply a proliferation of the vortex- 
skyrmions induced by disorder, because this finite p* is 
independent of the recursion equation of the disorder 
function R{z). This physically unaccepted p- value cer- 
tainly indicates that this fixed line is a unphysical one. 
In this way, within the model of the isotropic gradient 
energy, the quasi long-range order of the orbital orienta- 
tion, controlled by the nematic fixed point [s^ . is found 
to be stable against the vortex-skyrmions. 

To verify whether the above result is affected by the 
"orbital anisotropy" leading to the finite A and b, the 
same analysis as in the isotropic case will be applied to 
the full action ([55]) . Through lengthy but straightfor- 
ward calculations, we find that the one-loop recursion 
equations of the coefficients in ea. (l59p are given by 
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dT- 



dl 
d\2 

nr 

dl 

dp 

n 

dpi 

dl 
dp2 

dl 
dpi 

dl 

da 

H 
dui 

HT 

da2 

~dr 



f-\2- 


- e - JoR'^'\l)il ~ X3 - 


P + Pl+ P2 


+ 2(72) 


7 




(1)[A2(1 + A3 


+ P- Pi 


- P2- 20-2) 


- Pi], 




-JoR''^ 


(1)(3A3- A2 4 




P- Pi- h 


-2(72) 


- 2pi - p2 + 3a - 8(71 - (72 ), 




(l)[p(l + A3 4 


-P-Pi- 


- P2 - 2(72 ) -\ 








(l)[pi(3 + A3 


+ P- Pi 


- P2- 2(72 ) 


- 12a 




-JoR^^ 


(l)[/52(5 + A3 


+ P- Pi 


- P2- 2(72 ) 


+ (7l + 


10(72], 




(l)[^l(l + A3 


+ P- Pi 


— 07 — 2(79 ) 


— Pi — 


6(7], 


-JoR^^ 


(l)a(4 + A3 + 


P- /5i - 


P2 - 2(72 ), 






-JoR^^ 


(l)<7i(4 + A3- 


K /5 - /5i - 


- P2 - 2(72 ), 






-JoR^^ 


(1)<72(8 + A3- 


1- p - pi - 


- P2- 2(72 ), 







(71) 



where pi = 327r/9i, /52 = 8Trp2/d, pi — 327rpi, (7 = 
1287rV(l - 1/d), (71 = 87r2cri(l - 1/d), and (72 = 
167r^cr2/((i((i + 2)). This set of equations have the fol- 
lowing fixed points or lines : 

(i) A2 = A3 = /5 = pi = /52 = Pi = (7 = (71 = (72 = 0, 

(ii) A3 + p~ pi~- p2-2a2 = -8, a = ai = 0, 

(iii) A3 + P - Pl - P2 - 2(72 = -5, (7 = (7l = (72 = 0, 

(iv) A3 + p- pi- P2-2a2 = -4, (72 = 0. (72) 

Among them, the resulting fixed point values of A3 in the 
cases (ii) and (iii) are found not to satisfy the elastic sta- 
bility condition, eg. ([5^ . In fact, we obtain A3 = —104/93 
in case (ii) and —5/3 in case (iii), respectively. Thus, 
these cases are unphysical. Further, in the case (iv), we 
find that p and a are always negative using the elastic 
stability condition A3 -f 1 > 0. Thus, just as in the sim- 
ilar situation in the isotropic approximation, this case 
is also judged to be unphysical. In contrast, the linear 
stability of the nematic fixed point (i) is easily verified. 
Then, if working around this nematic fixed point, the 
analysis on the disorder function R{z) is the same as in 
Ref. j9] and will not be repeated here. Therefore, we 
reach again the conclusion that the only possible stable 
fixed point in T ^ limit is expressed as the nematic 
one with no vortex-skyrmions. This conclusion that the 
orbital anisotropy is irrelevant is quite reasonable, judg- 
ing from the fact that, even in the liquid crystal case [39|, 
the fixed point expression of the Frank energy (i.e., with 
A2 = A3 = 0) is the continuum version of the ferromag- 
netic Heisenberg spin model with no orbital anisotropy. 
Further, the above result that, at least at weak disor- 



der, all topological defects can be irrelevant at long dis- 
tances implies that the superfluid rigidity defined from 
the current-current correlation function remains finite, 
because pure Goldstone modes play no roles of destroy- 
ing superfiuidity. 



V. SUMMARY AND DISCUSSION 

In this paper, we have shown through calculation of 
free energy that, in the GL region outside the critical 
region, the disordered ABM state is lower in free en- 
ergy than other candidates of an equal-spin pairing state 
detected as the A-like phase in superfluid '^He in aero- 
gel. The local anisotropy characteristic of the aerogel 
structure plays essential roles in reaching this conclusion, 
because an anisotropy favors more anisotropic pairing 
states. If the scattering events are fully isotropic, a much 
stronger disorder is needed for another ESP state to be 
realized, although, then, Tc itself would be extremely low- 
ered or vanish. The absence or presence of the genuine 
long range superfluid order is not essential to a possi- 
ble change of pairing states: In a situation with a long 
range order destroyed over some temperature range due 
to the thermal fluctuation, the entropic term lowers the 
free energy of some locally ordered state. The vortex liq- 
uid regime in the superconducting vortex phase diagram 
[41I |4J] is its typical example. Similarly, even in the 
present case where a static randomness destroys a long 
range order, a free energy gain from the random-field 
term overcomes a cost of the elastic (gradient) energy 

iim. 
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P[bar] 


BW 


ABM P(planar) 


R(robust) 


24 


1.243 


1.245 


1.445 


1.630 


ZO 




1.192 


1.414 




34.4 


1.210 


1.155 


1.399 


1.578 


24 


1.267 


1.278 


1.473 


1.654 


28 


1.244 


1.227 


1.443 


1.621 


34.4 


1.233 


1.190 


1.428 


1.603 



TABLE I: j3n/l3o{T) value at T = T^o of each pairing state 
for l/(27rrcoT) = (upper half) and 0.065 (lower half). 



APPENDIX B 

In this Appendix, derivation of the free energy of the 
continuum version of the random Ising spin model 



5, 



1 



^(0W)2 + ±(V0(a))2 
^(a)^(6))2 



(Bl) 



In the present work, we have given one possible sce- 
nario of the globally isotropic disordered ABM state with 
a finite superfluid density [l^: The A- like phase is an 
elastic glass [41*1 and is in the ABM pairing state with su- 
perfluidity as well as in '^He in aerogels with an uniaxially 
stretched anisotropy over large scales [12,9] . An alterna- 
tive scenario will be the case in which disorder-induced 
topological defects including the vortex-skyrmions are 
pinned by the disorder itself at time scales seen in real 
experiments. In this case, a nonvanishing superfluid re- 
sponse may be observed. At present, however, it is un- 
clear whether these scenarios assuming globally isotropic 
samples are relevant to real systems or not. In our opin- 
ion, for further development of the present subject, it is 
necessary for experimentalists to clarify whether globally 
isotropic aerogel samples are truly available among those 
used in experiments. 
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APPENDIX A 



(sc) 



In Ref . ■ the strong couphng (SC) corrections, 50. 

and 5(3^^'^'^ , to the GL-quartic parameters (3j were exam- 
ined in details. Based on the results obtained there, we 
list here the estimated pressure dependence of /3n (N = 
BW, ABM, P, and R) in Table I. 

M 

The data in Table I show that, with increasing disor- 
der, the SC correction in the ABM case is weakened more 
remarkably compared with those of other pairing states, 
leading to a rapid narrowing of the temperature range of 
the A- like phase (see Fig. 3). Nevertheless, this effect is 
not quantitatively substantial at all and does not lead to 
replacement of the ABM state with other one including 
the planar or robust state. 



will be reviewed based on Ref. [28| . The analysis proceeds 
as follows. First, we divide (j) into its mean field, which 
is (0)mf in T < Tc and zero in T > Tc, and a fluctuation 
S(f>^°'\ Next, the fluctuation part in 5ising is assumed to 
be well approximated by the trial action 

= / E^aV(p) SU-P) SMp)- (B2) 

P n J) 



Then, when calculated according to eg. ([29)1 . the free en- 
ergy is well approximated by 



^ = 2^/^trln(g-(p)) + -^((5. 



(B3) 



with taking n limit at the end, where V is the vol- 
ume, and tr denotes here the trace over the replica in- 
dices. Finally, F is calculated in terms of the solution of 
the saddle-point equations: 



SF 



SF 



SGabip) 



SQaa{p) 



5F 



S{<f) MF 



0, 
0, 
= 0. 



(B4) 
(B5) 
(B6) 



The replica-symmetry breaking, which may not be negli- 
gible in the critical region [i^ , will not be considered for 
Qab- Then, we have 



where 



Gab = Gcip)Sa,b + /-i iGcip))'^ 



Gc{p) = 



A-l-p2 



(B7) 



(B8) 



The "mass" A of fluctuation Scj) will be determined 
through eq.(B4). The parameter /x is related to the glass 
order parameter, which is, by definition [23|, nonvanish- 
ing below Tc, and is determined by eq.(B5), while the 
average value ((/))mf follows from eq.(B6). Further, we 
focus only on the region outside the critical region in 
which 



I I «|A|<1 



(B9) 
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and assume u <C g <C 1. The latter relation is safely sat- 
isfied in superfluid ^He at weak static disorder. Then, /i 
simply becomes 2u((0)mf)^, and the free energy density 
/ is expressed in the form 



^((0>MF)' + f((0)MFr, 



(BIO) 



where [F{gc{p))] = JpFiGcip))- Below, A wih be de- 
noted as Ap (A/) in T > (T < T^). The free energy 
density /p in T > Tc simply becomes 



fp 



a: 



3/2 



fc (Bll) 



except a constant fc depending only on a momentum cut- 
off Pc, where Xp — tq + {3g — 2u)pc/ (Stt^). The first term 
of fp is the ordinary Gaussian fluctuation term leading 
to the singular behavior ~ (T — Tc)^^/^ of the specific 
heat. Under the condition (|B9|) . this Ap^^ term may be 
neglected together with the corresponding one in // given 
below. The free energy density // below Tc is /(A = A/) 
and takes the form 



,3/2 
A/ 

127r 



- ^mMFf + lmMFY+fc 



X2 

16.g 



1 



(ig-2u)[gc? + fc. (B12) 



where 



and 



A/ = 2(-ro-(35-2u)[g,]), 



2g' 



(B13) 



(B14) 



Although the present analysis takes account of fluctu- 
ation effects, the critical region is neglected. Neverthe- 
less, the expressions (IB11[) and (jB12[) ensure a continuous 
transition at Tc defined byAp = = — A//2. 



APPENDIX C 

In this Appendix, we explain why the planar state is 
not realized in the GL region. To do this, let us first ex- 
amine the gradient energy in the planar state. The sym- 
metry variable a^_i of the planar pair-field is expressed 
as 



V2 



(CI) 



lilj, and R^^i is the real rotation 



where 5jj = Sij 
matrix expressing the BW state. Below, this 1-vector 
expressing the local anistropy axis in the planar state 
will be represented in terms of the same Euler angles 
as those in the ABM state (see eg. ([52]) '). After substi- 
tuting eq.(Cl) into the gradient energy any term 
unaccompanied by ^Sjf. can be neglected in the present 
harmonic approximation, because the disorder term in 
i^sym depends only on the Euler angle 9i expressing 1. Al- 
though a close examination is necessary for a cross term 
like Rp,,k{y Rp.,m)5j 5j this term is found to depend 
only on the Euler angle 0; in the present harmonic ap- 
proximation. Then, the gradient energy related to the 
disorder term is simply 



|Amf, 

K2 



2KI + K2 



J dV(^(divl)2 + ((l. V)l)2 



-(1-curll)' 



(C2) 



By applying the present harmonic approximation to 
cq.(C2) again, the resulting harmonic elastic energy is 
found to be 2.4 times bigger than the corresponding one 
eg. (155)) for the ABM case. It means that the free en- 
ergy gain due to the quenched disorder in the planar 
state is smaller than eq. (j56|) in magnitude. Further, since 
/3p > /3a, the planar state cannot become stable through 
Tamp (see sec. III). Therefore, no possibility of realizing 
the planar pairing state due to the impurity disorder is 
expected anywhere in the phase diagram at least in GL 
theory. 
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